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Abstract. The origin of the spiral pattern of magnetic fields in disc galaxies is an open question. Comparison of the regular 
magnetic field orientation with the gaseous spiral arm pitch angles can tell us whether spiral shock compression is responsible 
for the magnetic spirals. We also wish to see whether the ridges of different components of the ISM show the large-scale, 
systematic shifts expected from density wave theory. We have developed a technique of isolating elongated structures in galactic 
^ . images, such as spiral arms, using anisotropic wavelets and apply this to maps of the CO, infrared and radio continuum emission 

' of the grand-design spiral galaxy M51. Systematic shifts between the ridges of CO, infrared and radio continuum emission that 

QO , are several kpc long are identified, as well as large variations in pitch angle along spiral arms, of a few tens of degrees. We 

■ find two types of arms of polarized radio emission: one has a ridge close to the ridge of CO, with similar pitch angles for the 

0^ ' CO and polarization spirals and the regular magnetic field; the other does not always coincide with the CO arm and its pitch 

, angle differs from the orientation of its regular magnetic field. The offsets between ridges of regular magnetic field, dense gas 

and warm dust are compatible with the sequence expected from spiral density wave triggered star formation, with a delay of a 
few tens of millions of years between gas entering the shock and the formation of giant molecular clouds and a similar interval 
. between the formation of the clouds and the emergence of young star clusters. At the position of the CO arms the orientation 

of the regular magnetic field is the same as the pitch angle of the spiral arm, but away from the gaseous arms the orientation of 
JL , the regular field varies significantly. Spiral shock compression can explain the generation of one type of arm of strong polarized 

; I ■ radio emission but a different mechanism is probably responsible for a second type of polarization arm. 

-4— » 
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' 1 . Introduction versal shape for galactic spirals". Deviations from the logarith- 

H , mic spiral may be caused by local disturbances or arise from 

^ Disc galaxies often display spiral patterns in their distribu- . i u i » » u .u i * t 

■ - - ° f J f f systematic global enects, such as the gravitational force or a 

tions of stars, gas, dust and magnetic fields. Spiral structure is u i *u * »u • i j v 

' ° ' ° ^ nearby galaxy or the presence of more than one spiral density 

present in both the (i/sfnl7Mf/on of the radio synchrotron inten- ^t:., . , inon, Tf»u i c u • 

^ ^ wave (Elmegreen et al. 1989 1. If the regular magnetic field spi- 

sity, and the orientation of the kpc-scale regular field, as mea- ^^j^ ^^.^^ ^^.^^^^ ^.^^ ^^^^^^^ ^pj^.^! ^^^^ ^^^^^ 

sured by the B-vectors of radio polarization (e.g. Beck 2005 1. • c u u u u i i 

F V & compression, the magnetic field orientation should be closely 

Comparing the location and orientation of the spiral patterns in wj^^ui i - i it a , ^-e 

^ ° F F related to the local spiral arm pitch angle. In order to quantify 

different components of the interstellar medium (ISM) can pro- , n »u j ff » • i r a u » a e 

^ V / F jjoyy yygjj fjjg difterent spirals are aligned a robust method of 

vide important information on the astrophysical connections • i i u i • a 

^ ^ ^ measuring local pitch angles is required. 

between interstellar gas, dust and magnetic fields in galaxies. 

„ , . ..... , M5 1 is probably perturbed by a recent encounter with its 

Usually the prominent arms in disc galaxies are treated as . , xT^,^^,^^r;,- 

, . , . , . . , , , . , , J- companion galaxy NGC 5195. Such interactions often result 

having a logarithmic spiral pattern, where the pitch angle of . . , ;. . , , , 

, . , • , ■ TT ■ j j AAj i in increased star formation rates, either localised or global, as 

the arm is constant along its length. However. Kennicutt J1981I I • . , ^ a ■ ^a ■ u ■ 



found that the logarithmic spiral form is no better a mathemat- 
ical description of real spirals than hyperbolic forms, and that 
"... small scale distortions preclude the possibility of any uni- 



tidal forces and spiral density waves compress the interstellar 
medium. The two spiral arms of M51 can be traced through 
more than 360° in azimuthal angle in numerous wavebands 
and several authors have investigated their structure. Elmegreen 
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ral modes. They concluded that M51 contains an inner and an 
outer system of spiral arms with a conjunction of resonances 
at 2'- 3' galacto-centric radius, triggered by the companion. 
Measurements of the pitch angles of the main dust lanes by 
Howard & Byrd ( 1990) and of the variation in star formation 
efficiency along the arms by Knapen et al. (1992i confirmed 
this double spiral structure. It is clear from a cursory visual in- 
spection that even the well defined — grand design — spiral 
arms of M5 1 have many variations in pitch angle along their 
length. 

Radial shifts between the locations of the crest of the spiral 
arms in different constituents are expected from density wave 
theory. Petit et al. C199 6 ) observed a 300-400 pc shift between 
Ha and UV arms. Rand & Kulkarni ( tP^ . Scoville et al. 
<2001> . and Tosaki et al. (2002i noted that the Ha arms are 
located outside the CO arms. Henry et al. ( I2003> found that 
the Paschen-a line emission is offset from the CO emission 
along part of the arms. All these shifts refer to the inner spiral 
system and are expected for arms triggered by density waves. 
Tilanus et al. J1988t found that thermal radio continuum and 
Ha arms are significantly shifted outwards from the dust lanes, 
whereas the non-thermal radio continuum arms occur just in- 
side the dust lanes. 

Most of the above analyses were made by overlaying the 
contours of one image on a grey-scale map of the other or by 
making azimuthal cuts through the data at different radii. Both 
methods, while showing the most significant local displace- 
ments in the position of the arms, make it difficult to see the ex- 
tent along the arm of any systematic shift. In the present paper 
we describe an objective method, using wavelet analysis, to de- 
termine the location of the spiral arm ridges, and measure their 
local pitch angles, along the entire length of an arm. We apply 
this method to images of the total and polarized radio contin- 
uum emission at A6 cm (Fletcher et al. in preparation), the dust 
emission at 15 /im (Sauvage et al. ll996> and the CO(l-O) fine 
emission (Heifer et al. 2003 1. 

In Section|2lwe describe the anisotropic wavelet transform 
and how we use this method to identify the location of a spi- 
ral arm and to measure its localised pitch angles. The obser- 
vational data we use are briefly discussed in Section |3l and the 
results of applying the anisotropic wavelet method to the data 
are given in Section |3 The astrophysical implications of the 
results are discussed in Section|5l 

2. The Method 

In image analysis, wavelet based methods allow the isolation of 
features, such as spiral arms, and the decomposition of a map 
into a hierarchy of structures on different scales. Wavelets are 
a tool for data analysis based on self-similar functions which 
are well localised both in the physical and frequency domains. 
Using one-dimensional isotropic wavelets Frick et al. (.2000.) 
identified systematic shifts between the magnetic and optical 
spiral arms in NGC 6946 and using two-dimensional isotropic 
wavelets Frick et al. (2001.) investigated the scale-by-scale cor- 
relations between maps of the same galaxy in different spectral 
ranges. In this paper a new method, using a two-dimensional 
anisotropic wavelet is presented. 



The wavelet transform can be considered as a general- 
ization of the Fourier transform. The classical Fourier trans- 
form is based on harmonic functions. The generalized Fourier 
transform allows using non-harmonic orthogonal functions as 
a basis (e.g. a Walch set of discrete, piecewise-constant func- 
tions). The short-time Fourier transform and the Gabor trans- 
form use oscillatory basis functions with local support. The 
wavelet transform also uses oscillatory functions, but in con- 
trast to the classical Fourier transform these functions decay 
rapidly toward infinity and all functions in the wavelet basis 
are self-similar (the main distinction between the wavelet trans- 
form and the Gabor transform). One-dimensional and isotropic 
multi-dimensional wavelet transforms are based on the space- 
scale decomposition of the signal (in other words, the family 
of wavelets has two parameters, governing the location and 
the size of the basis function). Using the continuous isotropic 
wavelet transform, a 2D image is decomposed into a 3D cube 
of wavelet coefficients. Cross-sections of the cube are slices 
which contain the image details at a fixed scale. As a result, the 
wavelet transform conserves the local properties of the image at 
all scales. If required, the original image can be reconstructed 
from the cube by summing over all scales (this procedure is 
called the inverse wavelet transformation). 

An anisotropic wavelet transform is the convolution of the 
image with a set of wavelets having different locations, sizes 
and orientations. Such a family of basis functions is generated 
by translations, dilations and rotations of the basic wavelet. 
Applying the two-dimensional anisotropic wavelet transform 
to an image generates a four-dimensional data set, which is 
a space-scale-orientation decomposition. Fixing the space and 
scale parameters — based on some objective criteria — en- 
ables one to track the orientation of an elongated structure 
(Antoine 1993). An extended description of the continuous 
wavelet transform can be found, for example, in Holschneider 
im^ i and Torresani ( fT995l . 



2.1. The Texan Hat function 

The ideal anisotropic wavelet for astronomical image process- 
ing should combine high angular sensitivity with a simple com- 
putational formula; the latter requirement is due to the time 
required to calculate the wavelet transform of a high resolu- 
tion image. In this work we use a specially designed ad hoc 
anisotropic wavelet. Its formula is very simple and the wavelet 
transformation can be calculated efficiently. We do not use the 
Cauchy wavelet (Antoine 1993 ), which provides the best angu- 
lar sensitivity, due to its more complicated formula. 

We introduce our wavelet by starting with the Mexican Hat 
function (MH), one of the most commonly used wavelets. In 
the ID case, the MH is described by i/'(x) - (l-x^)exp{-x^/2). 
The ID MH can be extended into 2D in one of two ways: 
(i) constructing an isotropic (axisymmetric) function by rota- 
tion of the ID MH, (ii) translation of the ID MH along an 
axis but restricting it within a (e.g. Gaussian-shaped) window 
in this direction, to obtain an anisotropic function. We call 
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the anisotropic wavelet constructed using the latter method the 
Texan Hat function (TH) and it can be represented by 



i//(x, y, a) 



1 



exp - 



(1) 



where a is the scale parameter by which the wavelet is dilated. 
In the form of Eq. Q the TH is sensitive to structures elongated 
along the x-axis. This means that if an image is convolved with 
the TH, extended objects parallel to the x-axis will be ampli- 
fied. Substitution of (x cos (p+y sin (p) for x and (y cos (p-x sin (f) 
for y in Eq. Q makes the TH sensitive to its orientation (p mea- 
sured from the x-axis, since this substitution is equivalent to a 
rotation of the axes by (/?. A complete set of basis functions of 
different sizes and orientations is obtained by dilation and rota- 
tion of the basic wavelet. Note that since the TH is a symmetric 
function with respect to rotation by 180° it is sensitive to the 
orientation of an elongated structure, not to its direction: 



>lf<p(x,y) = ^<p+m''ix,y). 



(2) 



Another, more formal, way to design the TH is as follows. 
The ID MH can be defined as the second derivative of the 
Gaussian function. Similarly, the TH is defined as the second 
partial derivative of the 2D Gaussian function: 



x^ +y^ 



(3) 



Equation can be parametrized for variation in its size a by 
substituting x with x/a and y with y/a (i.e. a dilation of the 
axes). This approach is useful for determining the properties 
of the TH in the Fourier domain. For example, one can easily 
estimate the Fourier transform of the TH by using the Fourier 
transform of the Gaussian function and the theorem of deriva- 
tives in the Fourier domain. Furthermore the definition of the 
wavelet as a derivative is not only valid for the 2D case, but 
also for higher dimensions. 

The anisotropy of the TH can be successfully exploited 
for an analysis of the orientations of galactic spiral arms. A 
brief overview of the method was first presented in Patrikeev et 



al. (I2005t . A detailed description of the technique and an illus- 
trative application to an artificial spiral image are given in the 
next section. 

2.2. Measurement of spiral arm position and pitch 
angles 

It is obvious that for an isotropic structure (e.g. a round-shaped 
object, or randomly distributed set of spots, etc.) an orientation 
cannot be defined. In this work we use the TH wavelet to trace 
the ridge of a spiral arm and to measure its local pitch angles. 
To illustrate how the TH can determine the position and pitch 
angle of a spiral arm, we use the artificial image presented in 
Fig. nja). The test image consists of: (1) two logarithmic spi- 
ral arms with Gaussian cross-sectional profiles, the first has a 
constant pitch angle of 20° and the second, broken, arm has a 
pitch angle of 15° at small radii and 25° at large radii; (2) an 
exponential disc that is about 70% of the arm intensity in the 



Fig. 1. Analysis of an artificial spiral image with the Texan Hat 
function, (a) Two arm logarithmic spiral with a central expo- 
nential disc and randomly distributed, slightly smoothed, point 
sources. The pitch angle of the first (unbroken) arm is 20°; the 
pitch angle of the second (broken) arm is 15° at small radii and 
25° at large radii, (b) Map of the maximum wavelet coefficients 
at a fixed scale that is close to the average arm width. 
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Fig. 2. Azimuthal variation of pitch angle determined from the 
anisotropic wavelet transform of the spiral arms of Fig^. Pitch 
angles are measured at fixed azimuthal increments along the 
ridges of the arms. The azimuthal angle is measured counter- 
clockwise from the top of the image. 



central region; (3) fifty bright point sources, slightly smoothed 
with a Gaussian kernel, distributed randomly over the image. 

The test image was convolved with a TH at a fixed scal- 
ing parameter a close to the average spiral width and a rota- 
tion parameter tp ranging from 0° to 180° in steps of 1°. From 
the convolution we obtained a stack of maps of wavelet co- 
efficients; each map represents a specific value of (p and has 
the same size as the test image. Figure ^b) shows the map of 
maximum wavelet coefficients, obtained by taking the largest 
value through the stack at each position. The central exponen- 
tial background is invisible because the wavelet coefficients of 
such a structure are very low. The point sources are converted to 
faint spots; their wavelet coefficients have no pronounced max- 
imum since a round source has no orientation. In contrast to 
the non-oriented structures, the two-fold spiral becomes more 
emphatic, with large wavelet coefficients. 
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We define the location of the ridge of a spiral arm as the 
position where the maximum wavelet coefficient is found. The 
local pitch angle is defined as the TH orientation correspond- 
ing to the maximum wavelet coefficient. Figure |2 shows the 
measured azimuthal variation of pitch angle at the ridge of the 
spirals (the spike-like deviations were suppressed). The first 
arm has a pitch angle close to the expected value of 20°. The 
pitch angle of the second arm changes from 15° to 25° (near 
the azimuth about -270°) remaining constant before and after 
the break. Small deviations of the pitch angle are caused by the 
influence of bright sources located in the vicinity of the spirals. 

The analysis of the test image demonstrates the capability 
of the proposed method to determine the pitch angle of the spi- 
ral and to resolve the difference of 5° despite the exponential 
background and small-scale bright sources. 

3. The Data 

M5 1 is a well observed galaxy and data in many spectral ranges 
are available. For our analysis we use four maps: CO(l - 0) line 
emission of molecular gas, total and linearly polarized radio 
continuum emission at A6.2 cm, and mid-infrared dust emission 
at ISfxm. 

The CO(l -0) line emission at A2.6 mm was observed with 
the BIMA interferometer and corrected for missing spacings 
with single dish observations from the NRAO 12m dish on Kitt 
Peak (Heifer et al. '20031. CO emission is the best available 
tracer of ffie molecular gas in M5 1 . 

It is often assumed that dust lanes are the best tracers of 
large scale spiral arm shocks, although recent Hubble Space 
Telescope images show a plethora of dust-absorption spurs and 
branches attached to ffie M51 spiral arms that clearly require 
a more subtle interpretation. The recent 850 fim map of M5 1 
at 15" resolution — twice the resolution of the data used in 
ffiis paper — by Meijerink et al. J2005> shows that most of 
the (cold) dust in M51 originates from an exponential disc, 
rather than from narrow spiral arms associated with the dust- 
absorption lanes seen in the optical bands. We note that at high 
resolution the strongest CO(l - 0) emission appears to closely 
follow the main spkal dust lanes in the inner r < 4 kpc (Aalto 
et al. [T999l . 

The /i6.2cm total and linearly polarized radio continuum 
emission was observed with the VLA ' and corrected for miss- 
ing short baselines using maps obtained with ffie Effelsberg ^ 
telescope (Fletcher et al. in preparation). The total emission 
is a mixture of synchrotron radiation of cosmic-ray electrons 
spiralling in interstellar magnetic fields and bremsstrahlung 
emission from thermal electrons. The relative fractions of syn- 
chrotron and thermal emission vary with position; at the lo- 
cation of H II regions in the spiral arms the emission is up to 
50% thermal but in the inter-arm regions this fraction falls to 
< 20%. Hence, the total radio intensity 16 is mostly a mea- 
sure of the total strength of the magnetic field and the cosmic 

' The VLA is operated by the NRAO. The NRAO is a facility of the 
National Science Foundation operated under co-operative agreement 
by Associated Universities, Inc. 

^ The Effelsberg telescope is operated by the Max-Planck-Institut 
fiir Radioastronomie on behalf of the Max-Planck-Gesellschaft. 



ray electron density. The polarized emission PI6 is purely non- 
thermal, and its intensity is a measure of ffie strength of the 
regularly oriented magnetic field. The polarization B-vectors 
were corrected for Faraday rotation using new /I3.5cm data 
(Fletcher et al. in preparation) and show the orientation of the 
regular magnetic field. The similarity of the polarized intensity 
maps at AA3.5, 6.2 cm (apart from generally reduced intensities 
at A3. 5 cm in line with the non-thermal spectral index) and ffie 
low Faraday rotation measures mean that Faraday depolariza- 
tion does not significantly effect the PI6 map. 

The 15 fim dust emission (denoted ISO throughout ffiis pa- 
per) was observed with ISOCAM (Sauvage et al. 1996|l and is 
a combination of thermal continuum emission from dust parti- 
cles and line emission from PAH molecules. 

All maps were smoothed to 8" resolution, equivalent to 
- 400 pc at the assumed distance of M51 (9.7 Mpc, Sandage & 
Tammann 1974 1. The maps were rectified to a face-on orienta- 
tion using an assumed inclination of 20° (where 0° is face-on) 
and a position angle of -10° for the orientation of the major 
axis (Tully 1974 1. Figure|3shows the face-on maps at 8" reso- 
lution. 

If the de-projection to face-on orientation uses incorrect pa- 
rameters and the spiral arms have a simple mathematical de- 
scription, such as logarithmic spirals, sine-wave like oscilla- 
tions will appear in log(radius)-azimuth plots of the ridges (see 
Fig. 2 of Kennicutt 1981 1. For r < 6 kpc we do not find any 
such oscillations in either arm, but for r > 6 kpc a series of 3 
peaks and troughs is present in the 16 and ISO of arm 1 (Fig.|6jl. 
The pattern may indicate a distortion of the disc in this region 
of M51, rather than intrinsic variability in the spiral pattern. 
Without detailed information on the velocity field we cannot 
investigate this possibility; however, whatever ffie cause of the 
oscillations, the relative positions of the 16 and ISO arms will 
be the same. 



4. Results 

4.1. The anisotropic wavelet transform 

Figure 0] shows the maximum wavelet coefficients for the four 
maps in Fig.|3l The scale parameter a of the wavelet was fixed 
to the approximate width of the spiral arms; tests showed that 
the results are insensitive to the aspect ratio of the wavelet in 
the range 0.5-2 times the chosen aspect ratio. The anisotropic 
wavelet transform clearly picks out the elongated arms in the 
images. In addition, smaller scale structures such as spurs 
roughly perpendicular to ffie arms are identified. In the PI6 map 
the low signal to noise ratio (relative to the other maps) and ffie 
patchiness of the diffuse emission results in a web of structures 
showing up in Fig. |3d), however extended arms of polarized 
emission are also evident. 



^ A note on terminology is required. In this paper we use the term 
"regular magnetic field" to refer to the magnetic field that gives rise 
to the observed polarized synchrotron emission. Polarized emission 
can be produced by an anisotropic random magnetic field as well as a 
coherent mean field; in this work we ignore the distinction except for 
a brief discussion in Sect.|5| 
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Fig. 3. The maps analysed in this paper, a) CO(l-O) emission (Heifer et al. l2003> . b Total radio emission at A6 cm (Fletcher et al. 
in preparation), c) Infrared emission at 15 fim (Sauvage et al. 1996 1, d) Linearly polarized A6 cm radio emission (Fletcher et al. in 
preparation), with B-vectors showing the orientation of the regular magnetic field (corrected for Faraday rotation). A grey-scale 
proportional to the square root of the intensity has been used to emphasise the spiral structures. Arms 1 and 2 are indicated on 
the CO map. The resolution is 8" in all maps and all have been rectified to a face-on orientation. The co-ordinate system is in 
seconds of arc relative to the galaxy centre (RA, Dec (2000):13h 27m 46s, +4-7° 27' 10")- 



The location of the wavelet coefficient maxima at a given 
radius are interpreted as the position of the spiral arm ridge. To 
ensure an even spacing of measured positions we select an ini- 
tial point that clearly lies on a spiral arm and then choose the 
location with the maximum wavelet coefficient that is a fixed 



increment in azimuth away from the initial position (2 degrees) 
within a arc of fixed opening angle. The anisotropic wavelet 
will pick out spurs and ridges that are connected to the spiral 
arms (Fig.0} but by selecting the maximum wavelet coeflicient 
within an arc we attempt to pick out the continuous structure of 
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position relative to centre [kpc] position relative to centre [kpc] 

Fig. 4. The maximum anisotropic-wavelet coefficients of the original maps shown in Fig.|3] corrected for inclination of the galaxy, 
a) CO(l-O) emission, b) total radio emission at A6 cm. c) infrared ISO 15 yum emission, d) polarized A6 cm radio emission. The 
grey-scale is proportional to the amplitude of the wavelet coefficient (white is high). At each location the wavelet coefficient was 
calculated for a fixed range of scales and 180 orientations, the highest coefficient is shown so the wavelet scale may be different 
between pixels. 



the spiral arm. The pitch angle of the arm is given by the orien- 
tation of the anisotropic wavelet that generates the maximum 
wavelet coefficient at each of these positions. 

Uncertainties in the measured pitch angle are caused by (i) 
the form of the real structures (they can be curved, asymmetric, 
having varying width etc.) and (ii) by instrumental noise in the 
map. In our approach structures on scales smaller than the scale 
of the arms - such as discrete bright clouds, spurs and breaks - 
can be considered as noise. To calculate the errors in our mea- 
sured pitch angles we used a Monte-Carlo technique. We define 
all structures on scales less than the approximate width of the 
spiral arms of 1 kpc as belonging to the noise. We carry out 
large-scale wavelet filtration of the image, keeping only scales 



a < I kpc and consider this to be a map of the noise. Then we 
reanalyse the original map overlaid with the noise map shifted 
(the noise map can be laterally shifted, rotated or reflected) by 
an arbitrary amount and calculate the pitch angles again. The 
process is repeated using several arbitrarily shifted noise maps 
to obtain the standard deviation for the pitch angle of the arms. 

4.2. The location of the spiral arms 

Figure |5] shows the location of the spiral arm ridges in M5 1 
derived from the anisotropic wavelet transform (Fig.E) of the 
total and polarized A6 cm, CO and ISO 15 //m maps. An alter- 
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Fig. 5. Location of the various spiral 
arm ridges in the plane of M51; dot- 
ted lines show the galacto-centric radii 
3, 6 and 9kpc. 16 and PI6 are the to- 
tal and polarized radio intensities at 
/16 cm, ISO the 15//m infrared emis- 
sion. The ridges are the positions where 
the anisotropic wavelet coefficient has 
a maximum, tracking along the spiral 
arms, as described in the text. 



native view of the same data — in log(r)-azimuth co-ordinates 
— is shown in Fig.|6l 

The small-amplitude 'wiggles' in the plotted ridge lines are 
due to the fixed number of positions where the ridge position 
was measured and do not translate into pitch angle fluctuations. 

There are three systematic offsets between the spiral arms 
traced in different wavebands; the most obvious offset is that of 
the polarized emission, but since this is the most complicated 
to discuss we will leave it to last. 

First, the CO spiral arms are consistently situated on the 
concave side of the ISO 15//m arms. The shift between the 
two ridges is typically ~ lOOpc, increasing to around 600 pc 
in the northern part of arm 2. This offset fits neatly with the 
chronology of the conversion of gas into stars in the large-scale 
shock model of spiral arms (Roberts 119691 Shu et al. 119721 
Tosa I1973> . The CO emission traces molecular gas clouds, 
whose formation is triggered by a shock on the upstream side 
of the spiral arm (the concave side within the co-rotation radius 
where the gas rotates faster than the spiral pattern). After dense 
cores in these clouds have collapsed, stars form downstream of 
the position where the clouds formed. The ISO 15 //m emission 
comprises line emission from PAH molecules and a continuum 
component from warm dust particles heated by Ly-a and UV 
radiation, the prime source of which are young star clusters that 
will have moved further downstream of the spiral arm shock. 
Figures|5]&|6lclearly show, for the first time that we are aware 
of, the systematic offset of dense gas spiral arms and infrared 
arms along regions several kpc long. 



Second, the 16 spiral arms are generally offset on the con- 
cave side of the ISO arms, by 200-300 pc. There is also a ten- 
dency for 16 to sit upstream of the CO arms inside ~ 6 kpc 
radius, by around 1 00-200 pc. This offset is harder to explain. 
Thermal radio emission will be enhanced by the same new star 
clusters that heat the dust and PAHs. Hence no offset between 
the 16 and ISO arms is expected from thermal effects. The most 
significant non-thermal process (on ~ lOOpc scales) associated 
with star formation is probably an increase in total magnetic 
field strength, caused by turbulent tangling of magnetic field 
lines. At our 8" resolution (~ 400 pc) this synchrotron emis- 
sion will be predominantly unpolarized, again tending to align 
the 16 and ISO ridges. The shift of the 16 arms upstream of the 
ISO arms can be due to the polarized synchrotron emission. 

The polarized ridge in the inner part of arm 2 is clearly lo- 
cated in the inter-arm region and may have no direct connection 
with the gaseous arms (see below), but along most of the inner 
part of arm 1, r < 6 kpc, the polarized ridge is close to the CO 
ridge, sometimes very closely aligned but often shifted about 
300 pc upstream. This displacement of the PI6 arm 1 indicates 
that the regular magnetic field is strongest upstream of the CO 
arm, and the PI6 ridge is shifted further upstream than the 16 
ridge. At radii of 6-8 kpc the PI6 ridges of arm 1 and arm 2 
move from the concave to the convex side of the other tracers. 
The shift was also seen in the /1/I18,20cm polarization maps 
of Horellou et al. J1992> . Beyond this radius, the two polariza- 
tion arms behave differently: along the outer part of arm 2 the 
polarization ridge is intertwined with the 16 and ISO, whereas 
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Fig. 6. Radius (on a logarithmic scale) versus azimuth plots of the spiral arm ridges. This shows the regions of systematic offsets 
in ridges quite well (compare to Fig.|5}. The azimuthal angle is measured counter-clockwise from the top of the face-on galaxy. 



along the outer part of arm 1 the PI6 sits consistently on the 
convex side of the 16 and ISO ridges. CO is absent at these 
large radii. The differences between the inner and outer spirals 
are difficult to explain, but may be related to the conjunction of 
an inner spiral density wave and an outer material spiral pat- 
tern identified by Elmegreen et al. ( 1 989j . overlapping in the 
range 6 kpc < r < 8 kpc. The largest shifts between PI6 and 
other tracers occur on the side of M5 1 closest to the companion 
galaxy, which may have caused the two outer spirals (Howard 
&Byrd[l990i. 

The offset of the PI6 ridge upstream of the CO ridge in 
arm 1 in the inner disc indicates that the spiral shock lies 
slightly upstream of the CO ridges. The position of the 16 ridge 
is then due to a balance between the polarized (non-thermal) 
emission pulling it upstream of the CO and the thermal and 
unpolarized emission pulling it towards the ISO ridge. 

Third, radio polarization arms are partly situated in inter- 
arm regions, most noticeably over ~ 3 kpc in the inner part of 
arm 2 and north of the middle part of arm 2 and in the outer 
part of arm 1, where they are displaced by around 700 pc from 
the other arms. Polarized intensity maps at XXh, 6 cm (Fletcher 
et al. in preparation) clearly show concentrations of emission in 
these inter-arm regions, particularly in the inner part of arm 2: 
we emphasise that these elongated polarized structures are not 
artefacts caused by the wavelet analysis. These sections of the 
PI6 ridges resemble the magnetic arms observed in NGC 6946 
(Beck & Hoernes il996j . However, the PI6 ridge that is asso- 



ciated with arm 2 near the galaxy's centre forms a continuous 
structure that is initially well aligned with the CO ridge then, 
as the radius increases, moves into the inter-arm region, rejoins 
arm 2 and then crosses inter-arm space again to link up with 
arm 1 . This interlacing of large scale magnetic spiral structures 
and gaseous arms is quite unlike the behaviour of the mag- 
netic arms seen in NGC 6946, but resembles that in NGC 2997 
(Han et al. 1999 1. Higher resolution and sensitivity may reveal 
similar structures in other galaxies. Its explanation will require 
more sophisticated modelling of the interplay between spiral 
density waves and magnetic fields than has been used to date 
and is beyond the scope of this paper. 

4.3. Variation of pitch angles along the spiral arms 

Figures and |8l show how the pitch angles of the spiral arms 
vary along their length. The uncertainties in pitch angles tend 
to get greater with increasing distance along the arms; the gen- 
erally fainter emission (compared to the inner disc) means that 
the maximum wavelet coefficients are smaller and the arms be- 
come less well defined (Fig.l^. First we will discuss the pitch 
angles for each spiral arm tracer in turn, then we will compare 
the pattern of pitch angles in different tracers. 

The two CO spiral arms can be traced through A(p » 360° in 
azimuth. In both CO arms the pitch angle in the inner quarter 
turn (A0 ^ 90°) is approximately constant, within the errors, 
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Fig. 7. Pitch angles of arm 1 (see Fig.|3lfor arm labelling), a) 
CO(l-O) emission, b) 16 A6cm total emission, c) ISO 15/vm 
emission, d) PI6 A6cm polarized emission. Every 6° in az- 
imuth (measured counter-clockwise from the top of the face-on 
galaxy) along an arm, the pitch angle of the wavelet with the 
highest amplitude is shown. Errors are calculated according to 
the Monte Carlo method described in Section BTTI The arms run 
from left (smallest radius) to right (largest radius) in all panels. 



but the pitch angles of the two arms in this region differ by 
~ 10°: in arm 1, pco ~ 25° and in arm 2, pco ~ 15°. Further 
along the arms the pitch angles follow different trends. Thus, 
despite the apparent similarity of the two CO arms, they are 
not symmetric under rotation. There is substantial variation of 
pco in both arms, with ranges 0° < pco ^ 30° in arm 1 and 
-20° < pco ^ 40° in arm 2. There is a tendency for low or 
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Fig. 8. As in Fig.0but for arm 2. 



negative pitch angles to occur toward the ends of both arms. 
This shift in pitch angles to low and even negative values oc- 
curs at ^ -100° in arm 1 and <p « +60° in arm 2. The tran- 
sition occurs around the place where the character of the CO 
arms - especially arm 2 - markedly changes in the original map 
(Fig-E)- The galacto-centric radii of these positions are around 
6 kpc, which is the approximate position of the co-rotation ra- 
dius (Elmegreen et al. ll992> . At these azimuths the pitch angles 
of the other tracers show similar variations. Furthermore in the 
azimuthal ranges 40° > > -100 in arm 1 and 200° >(f> > 60° 
in arm 2 the maxima and minima in pco follow each other with 
the same interval in 0; changes in pco occur symmetrically 
along the two arms, but with different magnitudes. 



10 



Patrikeev et al.: Gas, dust and magnetic spiral arms in M51 



The A6 cm total radio emission (16) and ISO 15 jum spirals 
are both quite well defined over the first half turn (A(p ^ 200°) 
with systematic trends in their pitch angles. The variation in 
pitch angles is stronger than in the CO arms and the uncer- 
tainty in the outer arm pitch angles is higher; this is because 
the arms are better defined in CO, whereas in 16 and ISO the 
arms are broader and contain bright patches, probably arising 
from regions of intense star formation. The A6 cm polarization 
arms (PI6) are generally the least well defined of all, although 
the first A(f> ^ 130° of arm 2 is very clear and shows a sharp and 
systematic variation in pitch angle from ppi ^ 30° to ppi ^ 0° 
and back up to ppi ^ 60°. Both PI6 arms show strong variations 
in pitch angle - a range of -20° < ppi < 60° in arm 2. 

In Section l4!2l we saw that the position of the spiral arm 
ridges is not the same in all tracers, but that there can be small 
but systematic shifts in the arm location. Now we consider how 
the pitch angles compare along arms. 

In general there is quite good agreement and pco ~ pi, 
particularly in arm 1 for (p > -120° where the pitch angles 
agree within errors for Acp ^ 100°. In arm 2 there are alternating 
regions of good and bad agreement in pitch angles. Overall, 
the similarity in pitch angles and locations where both tracers 
show changes to higher or lower pitches, the physical process 
producing the CO arms is likely the same as that producing the 
16 arms, as expected for density wave compression. 

Finally we compare pi and /?iso- In arm 2 the angles are 
for the most part equal, and where they are not the tendency 
(i.e. whether there is a decrease or an increase in pitch angle) 
is the same. Thus, despite the systematic small offset in the lo- 
cation of the 16 and ISO arms, their pitch angles are apparently 
modified for the same reasons. 

4.4. Orientation of the regular magnetic field and pitch 
angles of the gaseous and magnetic arms 

In this Section we will compare the orientation of the regular 
magnetic field lines with the pitch angles of the gaseous spiral 
structure (traced by CO emission) and with the pitch angles of 
the spiral structure in polarization (PI6). 

The intrinsic orientation of the regular magnetic field (/7b) 
was determined from the observed Stokes Q and U values 
by using observations at /I3.5cm and /I6.2cm to correct for 
Faraday rotation. The errors in magnetic field pitch angles were 
calculated using the noise in the observed maps only and do 
not include any systematic effects, such as uncertainty in the 
galaxy inclination (such systematic effects will equally effect 
the wavelet derived spiral arm pitch angles and so should not 
invalidate our comparisons). We only plot p^ where the polar- 
ized intensity at both wavelengths is greater than three times 
the noise level. In order to increase the signal to noise ratio and 
thus to see the azimuthal trends in ps better, we have used Q 
and U data smoothed to 15" (the anisotropic wavelet effectively 
smoothes the original 8" data in seeking the optimum position 
and orientation of the spiral arms). However, we have also cal- 
culated pb at 8" resolution and the main results — concerning 
the alignment or misalignment of pB with pco and ppi — do 
not change. 



Figure |5] shows how the orientation of the magnetic field 
lines varies with azimuthal angle at radii in the range 0.9 kpc < 
r < 6.5 kpc. The pitch angles of the two CO spiral arms are 
also indicated (see Sect. 14. 1> . These show the positions of the 
gaseous spiral arms and let us measure how well the orienta- 
tions of the magnetic field lines and spiral structure agree. 

The first thing to note is that the orientation of the mag- 
netic field lines is almost never zero. This means that the reg- 
ular magnetic field has a spiral shape — and therefore both 
a radial and an azimuthal component — virtually everywhere 
in the region shown. The orientation of the regular magnetic 
field is consistently less than 45° at the radii shown (pB > 45° 
only for small regions at r > 5.3 kpc); the azimuthal compo- 
nent of the magnetic field dominates the radial component for 
r < 6.5 kpc. 

Figure|5]also shows that the orientation of the regular mag- 
netic field lines vary by at least ~ 30° in each ring; this is by no 
means obvious in maps showing the B-vectors (e.g. Fig.|3ll), 
where the field orientation appears rather constant. The same 
situation may also exist in observations of other external galax- 
ies, whereby the field appears to be better aligned with the op- 
tical spiral and have a more consistent orientation than is really 
the case. Note that even for galaxies with a small inclination to 
the line of sight a slice at constant radius is quite elliptical and 
it is hard to judge "by sight" how the orientation of the field 
varies with azimuth. 

The other striking feature of Fig.|9]is the extremely good 
agreement between the CO spiral arm pitch angles and the 
magnetic field orientation at the position of the CO arms. 
This cannot be a coincidence as pB itself varies considerably 
at each radius and the agreement holds whether pco is low 
(e.g. Pco ~ 20° at r = 1.6 kpc) or high (e.g. pco ~ 40° at 
r = 5.0 kpc). Thus the orientation of the regular magnetic field 
and the gaseous spiral arms are tightly linked in M51. There 
is not a consistent trend in pB away from the spiral arms. The 
magnetic field orientation is sometimes higher and sometimes 
lower in the inter-arm region than at the position of the gaseous 
spiral arms. 

Berkhuijsen et al. J1997t found a broad agreement between 
optical arm pitch angles (Howard & Byrd 1990), averaged in 
~ 5 kpc^ sectors, and their large scale magnetic field model of 
M51. 

In Fig.[lO|we compare pB with ppi. In polarization arm 1 
Pb - Ppu the regular magnetic field orientation is parallel to 
the ridge of the PI6 arm. This indicates that the same physical 
effect is responsible for both the increase in polarized emission 
along the ridge and the orientation of the regular magnetic field 
there. This section of polarization arm 1 lies 1 00-200 pc up- 
stream of the corresponding arms in other tracers. In contrast, 
along polarization arm 2, pB < ppi, the magnetic field always 
has a lower pitch angle than the polarization arm, with a dif- 
ference in angle of 10-20°. This section of polarization arm 2 
lies about ~ 400 pc upstream of the other tracers of arm 2 and 
proceeds to cross over the gas and optical spirals at a radius of 
r ^ 5.5 kpc. This indicates that different forces are at work in 
producing polarization arms 1 and 2. 
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Fig. 9. Azimuthal variation of the orientation of the regular magnetic field pB = 1/2 arctan(2/J/), corrected for Faraday rotation, 
plotted as curves at different radii from the centre of 1VI51. Uncertainties are shown by upper and lower curves. The resolution 
of the observations is 15"Angles are only plotted where the polarized intensity at both A3 cm and A6 cm is above 3x the noise. 
Also shown are pitch angles of the CO arms (filled squares with error bars — most errors are slightly smaller than the size of the 
squares). The azimuthal angle is measured counter-clockwise from the top of the face-on galaxy and the gas flow relative to the 
spiral pattern is from left to right. 



5. Discussion 

New observations of the spiral galaxy M51, obtained with the 
VLA and Effelsberg radio telescopes, and published maps of 
the molecular (CO) and mid-infrared (15jt/m) emissions, were 
used to compare the spiral patterns present in the magnetic 
field, dense gas and dust distributions. 

1. We have shown that anisotropic wavelet functions are 
appropriate tools to locate the positions of spiral arms in images 
of galaxies and to measure their local pitch angles. 

2. A systematic offset between the spiral ridges of molecu- 
lar gas (CO) emission and of mid-infrared (ISO 15 /urn) emis- 
sion can be followed over many kpc along the arms of M51 
(Figs. |5l&|6j. The typical offset is 100 pc coiTesponding to a 
typical time delay of ~ 10^ yr for gas well inside the co-rotation 
radius, assuming a circular rotation velocity of 200 km s ' 
(Garcia-Burillo et al. '1993 1, pattern speed of 38kms 'kpc ' 
(Zimmer et al. .2004.) and arm inclination of 20°. The delay 
can be interpreted as the interval between the formation of 
Giant Molecular Clouds (GlVICs) and the appearance of newly 
formed stars. A delay of ~ 10^ yr between dense gas accumu- 
lation and star formation was also estimated by Tosaki et al. 
J2002> from an analysis of the location of '^CO and '^CO emis- 
sion fines. 



We do not see a shift in the relative positions of the ISO and 
CO arms, i.e. the CO arm does not move from the concave to 
the convex side of the ISO arm, near to the co-rotation radius of 
r ^ 6.2 kpc determined by Elmegreen et al. (1992). This may 
be due to the non-trivial dynamical effects of two, overlapping, 
spiral arm patterns in M51 (Elmegreen et al. ll989> . 

3. We also found a systematic offset between the spiral 
ridges of total radio continuum emission and those of mid- 
infrared and CO emission, indicating a shift between regions 
of strong thermal and non-thermal radio continuum emission. 
While the thermal component is expected to closely follow the 
distribution of star-forming regions as traced by mid-infrared 
emission, the non-thermal component is probably enhanced 
due to compression of magnetic fields in a spiral shock front, 
located upstream of the mid-infrared spiral arm. Cosmic ray 
electron density is expected to be higher on the downstream 
side of the arms, where supernova remnants should be more 
common, but on scales of a few hundred parsecs this may not 
have a significant effect on the synchrotron intensity; we are as- 
suming that the distribution of cosmic ray electrons is mainly 
determined by the magnetic field, rather than the distribution of 
their sources. The purely non-thermal, polarized radio ridges 
he the furthest upstream of all the data as expected in this 
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Fig. 10. As in Fig.|9l lines show the orientation of the regular 
magnetic field, but the filled circles with error bars represent 
the pitch angle of the ridges of A6 cm polarized emission, ppi. 
The polarized ridge running from azimuth ~ 100° at the bottom 
to ~ 60° at the top is the prominent inter-arm polarization ridge 
in Fig.|5l 



scenario. This is in agreement with the results of Tilanus et 
al. (1988 1 who found a general arm sequence of: non-thermal 
emission - (cold) dust lane - thermal emission - Ho-. 

Inside co-rotation our observations reveal the following se- 
quence from upstream to downstream: polarized radio emis- 
sion, total radio emission, CO, infrared emission, interpreted 
as a result of the sequence of components: shock-compressed 
magnetic field, molecular gas clouds, and UV heated (warm) 
dust. This can arise if the spiral shock is located slightly ahead 
of the CO ridge, a plausible consideration if molecule for- 
mation and/or cloud condensation are triggered by the spiral 
shock. The compression of magnetic field, resulting in stronger 
polarized radio emission, occurs at the shock, but the forma- 
tion of dense CO clouds is only completed a few hundred par- 
sec downstream. Following similar arguments to those used in 
point 1 above, we can estimate the time delay between gas en- 
tering the shock front and the accumulation of dense GMCs: a 
typical polarized ridge offset of 300 pc corresponds to an inter- 
val of ~ 10^ yr. The total time between gas entering the spiral 
shock and the emergence of newly formed stars will then be a 
few tens of millions of years. 

The above discussion only applies to a diffuse component 
of the CO in the inter-arm region, which can be compressed 



by a large scale shock at the spiral arm. Molecular gas that 
is already gravitationally bound in clouds will not respond to 
the shock. Extended, faint CO(l-O) emission in the inner inter- 
arm region of M5 1 is clearly seen in the combined single-dish 
and interferometer data of Heifer et al. (2003 (, suggesting that 
at least a part of the molecular gas is diffiise. However the 
higher resolution, interferometer only, CO(l-O) maps of Aalto 
et al. ( 1999 1 reveal small clumps of inter-arm molecular gas, so 
that some — difficult to estimate — fraction of the upstream 
gas is probably already confined in clouds. This means that the 
shock strength will be difficult to infer from gas density ratios 
in the inter-arm and arm regions. On the other hand, Aalto et al. 
could identify steep velocity gradients both in and, particularly 
clearly, between giant associations of molecular gas which are 
in broad agreement with theoretical expectations of spiral den- 
sity wave shocks. The detailed physics of the gas response to 
possible shocks in M51 and the effect this has on the spiral 
arm morphology requires more sophisticated modelling than 
we can attempt here, especially the apparent upstream shift of 
the shock in the magnetic field; there are, however, strong in- 
dications that large scale spiral shocks occur in the vicinity of 
the gas arms and that these compress diffuse gas. 

4. The maintenance of both radial and azimuthal compo- 
nents of the magnetic field, and hence a non-zero pitch angle, 
is an important prediction of dynamo theory (Beck et al. '1 9961 
Shukurov 2000 1 but is very difficult to explain if the magnetic 
field is purely passive and stretched by shear in the galactic 
diff'erential rotation, since after a small number of galactic ro- 
tations the field would be completely circular, i.e. it would have 
a pitch angle ps - 0. 

5. Although the B-vectors of polarized radio emission from 
M51 seem to have the same orientations as the CO arms 
(Fig. 13, they smoothly change by about +15° around any 
fixed radius, without a systematic trend between the arm and 
inter-arm regions (Fig.|9}. Smaller pitch angles in inter-arm re- 
gions are predicted by some dynamo models and have been 
observed e.g. in NGC 6946 (Rohde et al. 1999 l. In M51, how- 
ever, the regular magnetic field structure is probably dominated 
by streaming motions and compression in the spiral shock (see 
below). The field is brought into good alignment at the CO arms 
or vice versa (Fig.|5) and then relaxes into a diff'erent configu- 
ration in the inter-arm regions that depends on e.g. the amount 
of compression in the arms, the energy of turbulent flows, the 
orientation of dynamo-generated field modes etc. 

6. We found a close alignment of the CO pitch angles and 
the regular magnetic field orientation at the position of the arms 
(Fig.|9j. This could be due to shock amplification of the regular 
field component parallel to the shock. 

If the field is frozen into the gas then bJ'^ - eB\\ and B*"^ - 
B±, where e - p^"^ Ip is the compression ratio in the gas density 
p and the superscript {s) labels quantities at the shock front. 
The deflection of the magnetic field A0, where is the angle 
between the field and the shock front upstream of the shock 
front, is given by 



A6» = e-ef'^ = arctanf — 
B« 



arc tan 



eBn 
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^1-1 



(4) 



where the last line holds for Bj_/B\\ ^ 1. For e - 4 (max- 
imum degree of gas compression by an adiabatic shock) and 
B±/B\\ - 0.5 (i.e. the magnetic field is initially inclined at 30° 
to the shock front, a conservative estimate c.f. Fig.|9jl we obtain 
AO ^ 20°. Thus a strong spiral shock can be expected to align 
the regular magnetic field rather well with the spiral arm. 

If the shock is weak (e < 4) or the upstream field makes 
a larger angle with the shock front then an additional source 
of alignment is required. The conversion of an isotropic ran- 
dom magnetic field component upstream of the shock into 
an anisotropic random field in the shocked region (where the 
anisotropic random field will be perfectly aligned with the 
shock) provides such a mechanism. Strong polarized emis- 
sion arising from anisotropic random magnetic fields formed 
in shearing shocks has recently been identified in the barred 
galaxies NGC 1097 and NGC 1365 (Beck et al. UM5\ . 

7. Some of the polarized radio emission (and hence some of 
the regular magnetic fields) forms ridges which are interlaced 
with the gaseous arms: the polarized ridge of arm 2 clearly lies 
in between gas arms 1 and 2 south-east of the central region. 
The magnetic field is not oriented parallel to this polarization 
arm but is inclined to the axis of the arm by up to 20° (Fig.llO>. 
In contrast to polarization arm 1 , spiral shock compression does 
not seem to be the origin of these enhanced regular magnetic 
fields. The field orientation seems to be between that of the 
polarization arm and that of the nearby CO arm 2, like finding 
a compromise between two different forces. One force causes 
the polarization arm to cross CO arm 2 to join the outer arm 1, 
the second tries to align the field with the gas arm. 
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